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Abstract 

The statistical theory of the dipole flexoelectric (FE) polarization 
in liquid crystals is used to calculate the temperature dependence of 
order parameters, the elastic constants and the FE coefhcients. Two 
systems with polar wedge-shaped and banana-shaped molecules are 
investigated. In both cases the FE coefhcients are proportional to the 
dipole moment component parallel to the molecule symmetry axis. 
The origin of the FE effect and microscopic pictures of the distorted 
phases are discussed. 

Keywords: nematic liquid crystals, flexoelectric polarization, elas- 
tic constants, statistical theory. 

1 Introduction 

In an ideal nematic liquid crystal the molecules are on average aligned along 
one common direction ±A^ |jy. There is usually some deformation of the 
alignment and the vector depends on the position r* in the liquid crystal. 
The free-energy density due to the distortion of the vector is expressed in 
terms of the vector derivatives and the elastic constants Ki [2] 

/, = ^K,{V ■ Nf + ^K,[N ■ (V X iV)]2 + ^Ks[N x (V x N)f, (1) 

where Ki, K2, and are splay, twist, and bend elastic constants, respec- 
tively. 
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In a deformed uniaxial nematic liquid crystal, there should appear in 
many cases a spontaneous dielectric polarization invented by Meyer [3] 

P = eiiV(V-iV) + e3(iV- V)iV 

= eiN{V ■ N) -e^N X {V X N), (2) 

where ei and 63 are the splay and the bend fiexoelectric (FE) coefficients, re- 
spectively. The appearance of the spontaneous polarization in liquid crystals 
as a result of orientational deformations is called the fiexoelectric effect. 

The FE coefficients and the elastic constants are macroscopic parameters 
that are important in practical applications of liquid crystals. Thus, it is 
interesting to study the connection between microscopic and macroscopic 
parameters in order to fasilitate new material production. In this paper we 
will focus on the FE effect. 

A microscopic mechanism of the FE effect was proposed by Meyer [3], who 
showed that it is a steric effect due to the asymmetry of the molecular shape. 
But Prost and Marcerou [4J noticed that the polarization in a deformed liquid 
crystal can be also produced as a result of a gradient in the average density 
of quadrupole molecule moments. In 1976 Straley developed a microscopic 
theory [3] by extending the Onsager theory of orientational ordering in hard- 
rod gases. Other mean- field theory was given by Derzhanski and Petrov [B]. 
The density-functional theory was derived by Singh and Singh were the 
dipole and the quadrupole contributions to the FE coefficients were taken 
into account. Somoza and Tarazona |8] tried to include the effects of the 
relaxation in the angular distribution function of the deformed nematic. As 
far as computer simulations are concerned, we would like to mention the 
papers by Stelzer et al. [9] and by Billeter and Pelcovits ^lOj, where pear- 
shaped molecules were studied. 

In this paper, we present two systems with wedge-shaped and banana- 
shaped molecules. In our calculations, we use a simple potential energy 
of interactions which allows us to obtain the macroscpic parameters of the 
systems and to understand the FE effect origin in the deformed phases. 

2 Expressions for fiexoelectric coefficients 

The microscopic free energy for the system can be derived in the thermody- 
namic limit (A^ —7- 00, V —7- 00, N /V = const) from the Born-Bogoliubov- 
Green-Kirkwood-Yvon (BBGKY) hierarchy \TT\ or as the cluster expansion 
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for the uniaxial systems [T^]. The total free energy F consists of the entropy 
term and the interaction term, namely 

F = Pent + Fint, (3) 

where 

f3F,nt = J drdRG{r, i?){ln[G(r, R)k] - 1}, (4) 

I^Fint = -^JdndRidf2dR2Gin,Ri)G{f2,R2)fi2. (5) 

Here fu = exp(— /3$i2) — 1 is the Mayer function, $12 is the potential energy 
of interactions, dR = d(j)d9 sin Odif), 13 = 1/ (ksT), and A is related to the ideal 
gas properties. The normalization of the one-particle distribution function 
Gis 

/ dfdRG{r, R) = N. (6) 

The equilibrium distribution G minimizing the free energy of Eq. ([3]) satisfies 
the equation 

ln[G(fi, R,)A] - I df2dR2G{f2, i?2)/i2 = const. (7) 

In the case of the homogeneous uniaxial nematic phase, the distribution 
function depends only on the molecule orientation G{f, R) = Gq{R). We will 
further assume that G depends on a single argument Gq{R) = Go{n ■ N), 
where determines the phase orientation. This is satisfied exactly for the 
molecules with Coov symmetry but we expect that it is plausible for the 
banana-shaped molecules with G2V symmetry. The molecule electric dipole 
moment is defined as 

jl = Hil + /i2m + /i3^, (8) 

where the unit vectors {l^m^n) describe the molecule orientation. 

In the case of the wedge-shaped molecules, the FE coefficients have the 
form pJij 

d = J dudRidR2fi2Go{ni ■ N)GQ{n2 ■ N)u^nun2x, (9) 
63 = IJ^3 dudRidR2fi2Go{ni ■ N)GQ{n2 ■ N)u^nixn2x, (10) 
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where n determines the molecule C^ov symmetry axis. In the case of the 
banana-shaped molecules, the FE coefficients can be written in the approxi- 
mated form provided that I determines the molecule C2V symmetry axis and 
n determines the long molecule axis [14j 

ei = iJ,i J dudRidR2fi2Goini ■ N)GQ{n2 ■ N)uJun2x, (H) 
63 = fii J dudRidR2fi2Go{ni ■ N)G'Q{n2 ■ N)uJi^n2x- (12) 

In both cases, the FE coefficients are proportional to the dipole moment 
component parallel to the molecule symmetry axis. It results from the sym- 
metries of the interactions and of the Mayer function. The elastic constants 
of the uniaxial nematic phase are ||15j 

= ^ 1 dudR,dR2fi2G'oini ■ N)G',in2 ■ N)uln^xn2., (13) 
PK2 = ^I dudR,dR2fuG'oin, ■ iV)G[,(n2 ■ N)uln^xn2., (14) 
PK^ = \j dudRidR2fi2G'o{n, ■ N)G'^{n2 ■ N)uln,^n2.. (15) 



3 Results 

The calculations were performed for the square-well potential energy of the 
form 

{+00 for (m/ct) < 1, 
-e for 1 < (u/a) < 2, (16) 
for (u/a) > 2, 

where u is the distance between molecules, r 2 — ri = uA, and cr depends 
on the molecule orientations and on the vector A. For cr, we used simple 
expressions with three parameters, ctq defines the length scale, (Xi defines the 
nematic term, and a2 defines the FE term. In the case of the wedge-shaped 
molecules a has the form 



a 



(A • ni)2 + (A ■ n2f + fT2(A ■ Hi - A ■ ^2), 



(17) 



whereas for the banana-shaped molecules 



(T = (To + O"! 



{A ■ n^f + {A ■ n2f + (r2{A ■ k - A ■ k 



(18) 
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We used the density NVmoiiy = 0.1, ai = 2ao, a2 = ctq. The molecule 
volume Vmoi was estimated from the mutually excluded volume. The FE 
coefficients were expressed in fJ^i/cTQ, the elastic constants in e/ao, and the 
temperature in e/ks- The potential energy of the form $12(^/0") allows us to 
express the FE coefficients and the elastic constants as a finite series of the 
order parameters where the coefficients of the expansion can be calculated 
analitically. As a consequence it is easy to change the set of parameters 
and to study different systems. 

The temperature dependence of the order parameters and the elastic con- 
stants is similar to the known results [13], therefore the FE coefficients will 
be discussed only. 

3.1 Wedge-shaped molecules 

In the system of the wedge-shaped molecules, on decreasing the temperature 
we meet the ffist order transition from the isotropic to the uniaxial nematic 
phase at Tq = 0.7287. The temperature dependence of the FE coefficients is 
presented in Fig. [H The physical picture is shown in Fig. [2J On decreasing 
the temperature, the bend coefficient 63 is almost constant while the splay 
coefficient ei increases monotonically. The wedge-shaped molecules ffi to 
the splayed structure and there is an excess of molecules pointing to the 
splay origin. The induced polarization has the opposite direction. In the 
bend structure the FE polarization points to the bend centre because the 
molecules are slightly turned in the opposite direction. 

3.2 Banana-shaped molecules 

In the system of the banana-shaped molecules, on decreasing the temperature 
we meet the ffist order transition from the isotropic to the uniaxial nematic 
phase at Tc = 0.6611. The temperature dependence of the FE coefficients 
is presented in Fig. [3l The physical picture is shown in Fig. HI The banana- 
shaped molecules ffi to the bend structure and the FE polarization points to 
the bend centre. In the splayed structure the molecules are slightly turned 
and the average polarization is opposite to the splay centre position. 
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4 Conclusions 



In this paper, the statistical theory was used to study the dependence be- 
tween the microscopic and macroscopic parameters of nematic hquid crystals. 
The temperature dependence of the order parameters, the FE coefficients, 
and the elastic constants was obtained for the two systems. For the wedge- 
shaped molecules the bend coefficient 63 is almost constant whereas the splay 
coefficient ei is changing monotonically. For the banana-shaped molecules 
the splay coefficient ei is almost constant whereas the bend coefficient 63 is 
changing monotonically. In both cases the FE coefficients are proportional to 
the molecule dipole moment component, parallel to the molecule symmetry 
axis. This is the long axis of the wedge-shaped molecules and the short axis 
of the banana-shaped molecules. It results from the symmetries of the inter- 
actions and of the Mayer function. In other theories, where the interactions 
are not explicitly used or they are more complicated, the dependence on the 
all dipole moment components is present. The results are consistent with the 
microscopic pictures from the paper by Meyer but two additional situations 
were taken into account, the splayed phase with banana-shaped molecules 
and the bend phase with wedge-shaped molecules. 

At present stage, the qualitative comparison between the theory and the 
experiment is very difficult. The experimental data on the FE coefficients 
are still scarce and sometimes contradictory [16j. On the other hand, a 
nontrivial dependence on the details of the chemical structure was shown 
[17J. A small chemical modification of the molecule can generate significant 
change of the flexoelectricity. Thus, we would like to give only the estimates 
assuming typical values of the model parameters. If we assume the molecular 
length (To = Inm, the interaction energy e — O.leV, and the electric dipole 
moment fi = ID, then we can estimate the values of the FE coefficients 
h/(Tq = S.SpC/m, the elastic constants e/ao = IQpN, and the temperature 
e/fcs = 11607^ (it is three/four times larger then the typical isotropic-nematic 
transition temperature). 

The theory should be improved in order to describe exactly real systems. 
For the distance a, the general expansion proposed by Blum and Torruella 
|18) can be applied in order to describe molecular shapes better. The Meyer 
function /12 used in the low density limit can be replaced with the direct cor- 
relation function C12 from the Percus-Yevick approximation. In many cases, 
liquid crystals are mixtures of different substances or different conformers 
and it should be taken into account. Finally, it is possible to consider the 
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FE effect in tfie biaxial nematic pfiase altfiougfi tfiis pfiase is still very rare. 
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Figure 2: The homogeneous, the splay, and the bend phases of wedge-shaped 
molecules. The solid lines denote the local direction A^. The z axis is the 
symmetry axis of the homogeneous nematic phase. In the splay phase there 
is an excess of molecules pointing down to the splay origin. In the bend phase 
there is an excess of molecules pointing to the left whereas the bend center 
is on the right. 
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Figure 3: The temperature dependence of the flexoelectric coefficients 
expressed in Hi/ctq in the case of banana-shaped molecules. The temperature 
T is expressed in e/ks- The I — N transition is at = 0.6611. 
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Figure 4: The homogeneous, the splay, and the bend phases of banana- 
shaped molecules. The solid lines denote the local direction A^. The z axis 
is the symmetry axis of the homogeneous nematic phase. In the splay phase, 
molecules are slightly turned and the average polarization direction is op- 
posite to the splay origin. In the bend phase, the molecule dipole moments 
point to the bend center on the right. 
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